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ABSTRACT 

We test the X-ray emission predictions of galactic fountain models against XMM-Newton measure¬ 
ments of the emission from the Milky Way’s hot halo. These measurements are from 110 sight lines, 
spanning the full range of Galactic longitudes. We find that a magnetohydrodynamical simulation of 
a supernova-driven interstellar medium, ■which features a flow of hot gas from the disk to the halo, 
reproduces the temperature but significantly underpredicts the 0.5-2.0 keV surface brightness of the 
halo (by two orders of magnitude, if we compare the median predicted and observed values). This 
is true for versions of the model with and without an interstellar magnetic field. We consider differ¬ 
ent reasons for the discrepancy between the model predictions and the observations. We find taking 
into account overionization in cooled halo plasma, which could in principle boost the predicted X-ray 
emission, is unlikely in practice to bring the predictions in line with the observations. We also find 
that including thermal conduction, which would tend to increase the surface brightnesses of inter¬ 
faces between hot and cold gas, would not overcome the surface brightness shortfall. However, charge 
exchange emission from such interfaces, not included in the current model, may be significant. The 
faintness of the model may also be due to the lack of cosmic ray driving, meaning that the model may 
underestimate the amount of material transported from the disk to halo. In addition, an extended hot 
halo of accreted material may be important, by supplying hot electrons that could boost the emission 
of the material driven out from the disk. Additional model predictions are needed to test the relative 
importance of these processes in explaining the observed halo emission. 

Keywords: Galaxy: halo — ISM: structure — X-rays: diffuse background — X-rays: ISM 


1. INTRODUGTION 

X-ray observations show that the halo of our 
galaxy contains hot, diffuse plasma. This plasma 
is observed both in emission, as a component of 
the ~0.1-1 keV soft X-ray background ( S XRB; 
e.g.. iKuntz fc Sirowdenl 120001: lYosliino et al.l 120091 : 
iHenlev et al.ll2010t iHenlev fc Sheltonl 1201311 . and in ab¬ 
sorption, in high- resolution X-ray spec t ra of bright back¬ 
groun d sources ( iNicastro et al.l l2002t iR.asmussen et ahl 


2QQ3I; lMcKeinan_et_alJ~ 




2006 


2007 


Bregman fc Lloyd-Dayiesl l2007t lYao fc Wand _ 

Yao et al.l l2009t iHagihara et al.l 20101 iGupta et al.l 
201211 . Although the Milky Way’s hot halo is well 


studied observationally, the details of the origin of 
this hot plasma remain uncertain. Understanding the 
relatiye importance to the hot halo of supernoya (SN) 
driyen outflows from the disk and inflows from the 
intergalactic medium is a key part of understanding the 
functioning of the Galaxy and its interaction with its 
envi ronment. 

In IHenlev et al.l (|2010L hereafter HIO), we tested mod¬ 
els of the hot halo plasma against 26 XMM-Newton 
observations of the high-latitude SXRB, by comparing 
the observed temperatures and emission measures of the 
halo with the distributions predicted by different phys- 


® Gurrent address: Department of Astronomy, Haverford Col¬ 
lege, 370 Lancaster Avenue, Haverford, PA 19041, USA 


ical mode ls. HlO’s analysis favore d SN-driven galactic 
fountains (l.lonng fc Mac Lo^l2006l hereafter JM06) as a 
major, possibly dominant, source of the hot halo plasma 
observed in emission, although these fountain models 
tended to overpredict the halo temperature. Additional 
support for the heating of the halo by disk SNe comes 
from the observation that some models of the halo’s 
global gas distribution (constrained by various obser¬ 
vational data) i mply that the hot halo may be convec- 
tively unstable (|Henlev fc Sheltonl 1201411 . However, HIO 
were unable to rule out the possibility that an extended 
halo of accreted m aterial also contributes to the emission 
((Grain et al.ll201(Tll . 

In this paper, we further examine the X-ray pre¬ 
dictions of galactic fountain models, in the light of 
two developments since HIO. First, we use a much 
larger set of measuremen ts of the Galactic halo emission. 
IHenlev fc Sheltonl ()201.1l hereafter HS13) measured the 
halo X-ray emission on 110 high-latitude XMM-Newton 
sight lines, an approximately fourfold increase over HIO. 
This is the largest set of measurements of the halo X- 
ray emission with GGD-resolution spectra assembled to 
date. Furthermore, these observations span the full range 
of Galactic longitudes, whereas HlO’s data set was re¬ 
stricted to I = 120°-240°. See Section[2]for a description 
of the observational data. (Note that HS13 discussed the 
energetics of galactic outflows versus extragalactic accre¬ 
tion as sources of the observed X-ray emission, but were 
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unable to distinguish between these two scenarios: both 
SNe and infall provide more than enough energy to power 
the observed emission, and either process could plausibly 
explain the observed variation of the surface brightness 
on the sky.) 

Second, it has been discovered that the JM06 fountain 
model included an unphysical inflow of hot gas from the 
vertical boundaries that adversely affected the model’s 
X-ray predictions (see Section El). We therefore exam¬ 
ine a new model of the SN-driven interstella r medium 
(ISM ), in which there is no such hot inflow (|Hill et al.l 
120121 hereafter H12). In addition, these newer simula¬ 
tions include results obtained with a non-zero interstellar 
magnetic field. See Section El for a description of these 
ISM models. 

Section m describes how we obtained the X-ray predic¬ 
tions from the ISM models. We present the results of 
the comparison of the model predictions with the obser¬ 
vations in Section El We discuss the results in Section El 
and conclude with a summary in Section [7l 


2. OBSERVATIONS 

We use HSlS’s measurements of the Galactic 
halo X-ray emission. They measured this emis¬ 
sion on 110 high-latitude XMM-Newton sight lines, 
selected from an all-sky XMM-New ton survey of 
the SXRB (jHenlev fc Sheltonl 1201211. HS13 ap¬ 
plied various filters to iHenlev fc Sheltoiil ’s (|2012h ob¬ 
servations in order to minimize the contamination 
from charge exchange (CX) emission from within 
the solar syste m, a time-var i able c ontaminant of 
SXRB spectra (ICravens et al.l 120011: IWargelin et ahl 


2004 ISnowden et ahl 2004: IKoutroumpa et al.l 2007 


Fuiimoto et al. 200'il : IHenlev fc Sheltonll2008l : lEzoe et al l 

20101 : iCarter et aLl l201lll ' In addition, HS13 ex¬ 

cluded certain features from their sample (the Scorpius- 
Centaurus superbubble, the Eridanus Enhancement, and 
the Magellanic Clouds). The HS13 data set contains 
measurements for ^4 times as many sight lines as HlO’s 
data set, spanning the full range of Galactic longitudes. 
Here we give a brief overview of HS13’s spectral mod¬ 
eling method and their halo results; see HS13 for more 
details, and for a comparison of their results with those 
from other recent studies of the SXRB. 

HS13 analyzed the SXRB spectrum for each sight line 
with a standard SXRB model, with components repre¬ 
senting the foreground, Galactic halo, and extragalactic 
background emission. The foreground emission was con¬ 
strained using_shndowmgdata from the ROSAT All-Sky 
Survey (ISnowden et al.ll2000ll . For all but one sight line, 
HS13 used a single-temperature (IT) collisional ioniza¬ 
tion equilibrium (CIE) plasma model to model the halo 
emission, obtaining the X-ray temperature and emission 
measure for the halo on each sight line. For the remain¬ 
ing sight line, HS13 added another, hotter component 
(with temperature T ~ 11 x 10® K), in order to model 
excess emission in the observed spectrum around ~1 keV 
(see their Section 3.1.2). For that sight line, we use the 
results for the cooler component (T ~ 2 x 10® K). 

HS13 detected emission from ^(2-3) x 10® K halo 
plasma on 87 out of 110 sight lines (79%), with a me¬ 
dian temperature of 2.2 x 10® kE| and a typical in- 


For some sight lines, HS13 were unable to constrain the halo 


trinsic 0.5-2.0 keV surface brightness of (1.1-2.3) x 
10“^^ erg cm“^ s“^ deg” . On the remaining 23 sight 
lines, HS13 give upper limits for the halo surface bright- 
ness. _ 

IHenlev et ahl (1201411 compared a subset of HS13’s re¬ 
sults with a measurement of the Galactic halo emission 
from an XMM-Newton observation of a compact shad¬ 
owing cloud, G225.60—66.40. The good agreement be¬ 
tween their measurement and th at from the nearest HS13 
sight line led iHenlev et ahl (|2014ll to conclude that HS13’s 
measurements are not subject to systematic errors, and 
can confidently be used to test models of the halo emis¬ 
sion. 


3. GALACTIC FOUNTAIN MODELS 

The JM06 and H12 SN-driven ISM simulations were 
carried out using FlashEI a parallelized Eulerian hydro- 
dynamical code with adaptive mesh refinement (AMR). 
In each case, the model domain was a tall thin box ex¬ 
tending to z = izinax, with peHodic boundary conditions 
on the vertical sides, and zero-gradient boundary condi¬ 
tions on the upper and lower surfaces. The model domain 
was initialized with gas in hydrostatic equilibrium. This 
gas was then heated and stirred stochastically by Type la 
and core collapse SN explosions, each with a frequency 
and z distribution appropriate for the Milky Way in the 
vicinity of the Sun. Each SN injected 10®^ erg of energy 
into a small region of the grid. 60% of the core col¬ 
lapse SNe occurred in clusters of seven to 40 explosions, 
while the remaining SNe occurred in isolation. The gas in 
the model domain was also subject to radiative cooling, 
and to diffuse heating representing photoelectric heat¬ 
ing of dust grains. The simulations were run at least 
long enough to eradicate the initial conditions. From 
our point of view, the most important feature of these 
ISM models is that the SN heating drives a fountain of 
hot (>10® K) X-ray-emissive gas into the halo. For more 
details of the models, see JM06 and H12. 

HIO tested the JM06 ISM model, which was carried out 
in a 0.5 x 0.5 x 10 kpc^ model domain, with Zmax = 5 kpc. 
HIO found that the X-ray emission measures predicted 
by this model were in good agreement with their obser¬ 
vations, leading them to conclude that galactic fountains 
are a major, possibly dominant, contributor to the hot X- 
ray emission in the XMM-Newton band (as noted in the 
Introduction). However, the JM06 model overpredicted 
the observed X-ray temperature. 

It has subsequently been discovered that the bound¬ 
ary conditions at the upper and lower boundaries of the 
JM06 model domain led to an un physical inflow of hot , 
high-pressure gas into the domain (|Mac Low et al.ir2012l : 
iJoiing et al.ll201^ H12). Early in the simulation, while 
the initial conditions were still being eradicated, hot gas 
from SNe moved upward through and eventually off the 
domain, causing the ghost cells just off the domain to be 
set to a high-temperature, high-pressure state. Subse¬ 
quent radiative cooling caused the halo pressure to drop, 
causing material to be drawn into the model domain. 
The state of this inflowing material was determined by 


temperature. In such cases, they fixed the temperature at 2.1 x 
10® K. 

^ Developed at the University of Chicago Center for Astrophys- 
ical Thermonuclear Flashes; http://flash.uchicago.edu/web/ 
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the state of the ghost cells, leading to a hot, high-pressure 
inflow. This inflow adversely affected the X-ray predic¬ 
tions derived from the JM06 model. 

The newer HI 2 model used outflow-only boundary con¬ 
ditions, and was carried out in a much larger domain— 
1 X 1 X 40 kpc^, with ^max = 20 kpc—and so does not 
suffer from the unphysical inflow problem of the JM06 
model. There are a few additional differences from the 
JM06 model. First, the Type la (core collapse) SN 
rate was slightly higher (lower) than that used in JM06, 
though otherwise the SN heating was the same. Second, 
the grid initialization was slightly different—JM06 ini¬ 
tialized their entire domain with gas at 10^ K, whereas 
H12 initialized their domain with gas at 1.15 x 10'^ and 
1.15 X 10® below and above \z\ ~ 1 kpc, respectively (the 
pressure was continuous across the interface). Also, H12 
employed a higher gas surface mass density than JM06: 
13.2 versus 7.5 Mq pc“^. Finally, H12 ran versions of 
their model that included a magnetic field—here, we ex¬ 
amine versions with (model bx50) and without (model 
bxO) a magnetic field. In model bx50, the magnetic field 
was initially horizontal and uniform in the xy plane, with 
a magnitude of 6.5 yG, but decreased with height such 
that the ratio of magnetic and gas pressures was con¬ 
stant. Note that the radiative cooling and diffuse heating 
were incorrectly applied in the original H12 simulations, 
but corrected models were described in their erratum. 
Here, we use results from the corrected simulations. 

It should be noted that HS13 found that the halo emis¬ 
sion measure tends to increase toward the inner Galaxy 
{I = 0°). However, because the models’ domains are tall 
thin boxes, we are unable to determine how the model 
predictions would vary with Galactic longitude or lati¬ 
tude. Instead, we test how well the models can reproduce 
the overall distributions of observed halo temperatures 
and surface brightnesses. 

4. CHARACTERIZING THE FOUNTAIN MODEL X-RAY 
EMISSION 

As in HIO, for a given model epoch we calculated halo 
X-ray spectra for 242 vertical sight lines, looking up¬ 
ward and downward from the Galactic midplane. The 
vantage points for these sight lines were arranged in an 
11x11 grid in the Gala ctic midplane, with grid sp acings 
of ss98 pc. We used the lRavmond fc SmithTdlQTTli spec¬ 
tral code (updated by J. C. Raymond & B. W. Smith, 
1993, private communication with R. J. Edgar) to cal¬ 
culate the X-ray spectra, assuming that the plasma is in 
CIE and is optically thin. We excluded material within 
100 pc of the midplane from the emission calculations, 
as such material is not in the halo. Note that the SXRB 
model used in the observational analysis (HS13) included 
a foreground component (in addition to the halo compo¬ 
nent, the results for which we use here; Section [2|). This 
foreground component accounted for the observed emis¬ 
sion from within ^100 pc of the midplane. 

The true halo emission is likely from plasma with 
a range of temperatures; in the observational analysis, 
this emission was characterized with a IT plasma model 
(HS13). Similarly, the emission predicted by the fountain 
models that we are examining here is from plasma with a 
range of temperatures. Therefore, to ensure a like-with- 
like comparison of the models with the observations, we 
characterized the predicted X-ray emission by creating 


synthetic XMM-Newton observations of the SXRB, and 
then analyzing the resulting spectra with the same SXRB 
model used in the observational analysis. This method 
is described in full in HIO; here we give an overview, and 
point out the differences from HIO. 

For each model sight line, we combined the predicted 
halo X-ray emission with models for the foreground emis¬ 
sion, the extragalactic background, the instrumental flu¬ 
orescence lines, and residual soft proton contamination 
(HS13). We folded the resulting spectrum through the 
XMM-Newton response function and added Poissonian 
noise corresponding to a typical field of view and expo¬ 
sure time from HS13. The simulations also took into ac¬ 
count the XMM-Newton quiescent particle background. 
We simulated a MOSl and a MOS2 spectrum for each 
model sight line, each of which we grouped such that 
there were at least 25 counts per bin (as in HS13). We 
then fitted the grouped spectra with the input SXRB 
model, but with the halo component replaced with a 
IT CIE plasma model. We use the resulting best-fit IT 
halo model to calculate the intrinsic 0.5-2.0 keV surface 
brightness, S'o. 5 - 2 . 0 - Thus, for each model sight line, we 
obtained an X-ray temperature and surface brightness 
which characterize the predicted X-ray spectrum, and 
which can be compared with the observed temperatures 
and surface brightnesses. Note that the model surface 
brightnesses obtained in this way were typically ^20- 
40% lower than those obtained directly from the model 
spectra. This is likely because the IT model used in the 
fitting cannot always accurately capture the entire model 
spectrum, calculated from a multitemperature plasma. 

We used the same foreground model as in HIO, but a 
different model for the extragalactic background. HI O 
used a single unbroken power-law (|Chen et al.l Il997[) , 
whereas we us ed the model from HS13: a double bro¬ 
ken power-law (jSmith et al.ll200^ rescaled to match the 
expected surface brightness of the so urces that fell below 
the source removal flux th reshold (jMoretti et al.l 120031 
iHickox fc Markevitchl200^ . This alteration in the extra- 
galactic model resulted in a change in the typical level of 
soft proton contamination; we adjusted our input model 
accordingly. The source removal flux threshold used in 
HS13 is lower than that used by HIO, resulting in more of 
the XMM-Newton field of view being excluded in the ob¬ 
servational analysis. We therefore lowered the assumed 
field of view for the simulated observations from 480 to 
410 arcmin^. However, we kept the assumed exposure 
time at 15 ks. 

In order to ensure that the simulated spectra had ad¬ 
equate signal-to-noise to constrain the IT halo model, 
we rescaled the input halo spectra such that they had 
a specified surface brightness. We undid this rescal¬ 
ing at the end, by dividing the output emission mea¬ 
sure by the same factor that was used to multiply 
the input spectrum. In HIO, we rescaled the spec¬ 
tra to give a 0.4-2.0 keV surface brightness of 2.06 x 
10“^^ erg cm“^ s“^ deg~^. Here, we found that the 
X-ray temperatures resulting from this procedure may 
depend weakly on the assumed surface brightness used 
to rescale the spectra. We therefore used three dif¬ 
ferent values to rescale the spectra: 0.5-2.0 keV sur¬ 
face brightnesses of 1.14 x 10“^^, 1.54 x 10“^^, and 
2.34 X 10“^^ erg cm“^ s“^ deg~^ (these are the quar- 
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Table 1 

Observed and Predicted Halo Temperatures and Surface Brightnesses, and Predicted Hot Gas Properties 










Properties of hot (T > 10^ i 

K) gas 



Obs. or Model 

Time 

Temperature 

*50.5-2.0^ 


("e) 


[KY'^ 


L 



(Myr) 


(10® K) 



(10 

cm“®) 

(10 

cm“®) 


(kpc) 


(1) 

(2) 


(3) 


(4) 


(5) 


(6) 


(7) 

1 

Observations 


2.22 

(2.01,2.64)® 

1.07 

(0.50,1.53)° 







2 

H12 bxO 

85 

1.72 

(1.19,2.38) 

0.016 

(0.002,0.062) 

0.43 

(0.19,0.79) 

0.55 

(0.28,1.02) 

0.21 

(0.10,0.36) 

3 

H12 bxO 

135 

1.67 

(1.10,2.34) 

0.006 

(0.000,0.031) 

0.35 

(0.15,0.66) 

0.45 

(0.17,0.85) 

0.16 

(0.04,0.34) 

4 

H12 bxO 

185 

1.98 

(1.62,2.50) 

0.020 

(0.011,0.068) 

0.11 

(0.09,0.13) 

0.15 

(0.12,0.25) 

3.22 

(2.91,3.50) 

5 

H12 bxO 

235 

1.74 

(1.33,2.73) 

0.012 

(0.003,0.035) 

0.042 

(0.036,0.050) 

0.078 

(0.046,0.123) 

7.98 

(7.05,8.94) 

6 

H12 bxO 

285 

2.22 

(1.77,2.59) 

0.023 

(0.010,0.058) 

0.055 

(0.043,0.070) 

0.098 

(0.064,0.150) 

7.09 

(5.32,9.98) 

7 

H12 bxO 

335 

1.88 

(1.68,2.17) 

0.010 

(0.006,0.040) 

0.037 

(0.035,0.044) 

0.059 

(0.040,0.112) 

9.76 

(8.29,11.17) 

8 

H12 bx50 

85 

1.28 

(0.93q.96) 

0.001 

(0.000,0.015) 

0.073 

(0.047,0.160) 

0.10 

(0.05,0.36) 

0.46 

(0.19,0.93) 

9 

H12 bx50 

135 

1.35 

(1.01,2.03) 

0.004 

(0.001,0.035) 

0.13 

(0.10,0.20) 

0.17 

(0.11,0.40) 

0.63 

(0.43,0.89) 

10 

H12 bx50 

185 

1.52 

(1.02,2.17) 

0.005 

(0.001,0.029) 

0.098 

(0.040,0.140) 

0.15 

(0.09,0.29) 

1.11 

(0.70,1.59) 

11 

H12 bx50 

235 

1.86 

(1.35,2.78) 

0.015 

(0.002,0.048) 

0.037 

(0.025,0.056) 

0.11 

(0.04,0.18) 

4.14 

(3.64,4.88) 

12 

H12 bx50 

285 

1.92 

(1.61,2.79) 

0.012 

(0.003,0.041) 

0.021 

(0.018,0.027) 

0.066 

(0.041,0.109) 

9.36 

(9.07,9.66) 

13 

H12 bx50 

335 

1.93 

(1.60,2.63) 

0.011 

(0.002,0.046) 

0.017 

(0.014,0.022) 

0.060 

(0.026,0.111) 

11.06 

(9.66,12.03) 


Note. — For each quantity, we have tabulated the median value, followed by the lower and upper quartiles in parentheses. For the model predictions, these 
quartiles were calculated from the sets of values obtained from the 242 model sight lines that we examined at each model epoch. Columns 3 and 4 contain halo 
temperatures and surface brightnesses, respectively. Columns 5—7 contain the mean electron densities, the r.m.s. electron densities, and the path lengths of the 
hot gas along the model sight lines, respectively. 

^ 0.5—2.0 keV surface brightness in 10“^^ erg cm~^ s~^ deg~^. 

^ Including only sight lines on which the temperature was free to vary. 

^ Latitude-corrected value (assuming a plane parallel halo geometry), including non-detections at their best-fit values. 


tiles for sight lines on which ~(2-3) x 10® K emission is 
detected; HS13, Table 2). 

We subjected the halo and extragalactic components 
of the model to interstellar absorption. The assumed 
column density does not strongly affect the results, but 
here too we decided to use the quartiles from HS13: 
1.26 X 10^0, 1.63 X 10^0, and 2.12 x 10^° cm-^ (cf. 
1.7 X 10^® cm“^ in HIO). Thus, each model sight line 
was characterized a total of nine times. For the compar¬ 
ison with the observations, we first combined the results 
obtained with the different rescaling surface brightnesses 
and column densities. 

5. RESULTS 

Figure [T] compares the X-ray predictions of the H12 
fountain model with HS13’s halo observations in the 
temperature-surface-brightness plane. We show predic¬ 
tions from several different epochs of the models (a) with¬ 
out and (b) with a magnetic field (models bxO and bx50, 
respectively). For comparison with HIO, Figure[Tfa) also 
shows predictions from a single epoch of the JM06 model 
(the latest epoch examined by HIO). 

Figure [2] compares the predicted temperature and sur¬ 
face brightness distributions of the H12 models with the 
observed distributions. Again, we show results from 
several different epochs of models bxO and bx50. The 
observed and predicted halo temperatures and surface 
brightnesses are summarized in columns 3 and 4 of Ta¬ 
ble [TJ respectively. The observations are in row 1, while 
the predictions from models bxO and bx50 are in rows 
2-7 and 8-13, respectively. Tabled] also summarizes the 
properties of the hot (T > 10® K) gas from each epoch 
of the models—columns 5-7 contain the mean electron 
densities, (rie), the r.m.s. electron densities, and 

the path lengths, L, of this gas along the model sight 
lines, respectively. Note that the r.m.s. electron density 
is more useful than the mean electron density for inter¬ 
preting the X-ray emission predictions, since the emission 
measure £ = {n4)L. 


The predictions from model bxO appear to undergo a 
slight oscillation in the temperature-surface-brightness 
plane. The median predicted X-ray temperature and sur¬ 
face brightness oscillate with peak-to-peak amplitudes of 
~0.3 X 10® K and ~0.3 dex, respectively, with a period 
of ~100 Myr. This oscillation may be related to the 
“bouncing” of the halo material reported by HI 2. 

At the earliest epochs of the bxO model shown here, 
the r.m.s. density of the hot gas is relatively high (~5 x 
10 “'^ cm“®) and the path length through this gas is rel¬ 
atively short (~0.2 kpc). At later epochs, the density is 
lower (~1 X 10“^ cm“®) but the path length much longer 
(up to ~10 kpc). However, these changes are such that 
£ remains the same within a factor of ~2 (~5 x 10“® and 
^1 X 10“^ cm“® pc, respectively, using the above-quoted 
densities and path lengths). (As an aside, we note that 
the median r.m.s. density and path length from the sight 
lines through the JM06 model are 9 x 10“^ cm“® and 
4 kpc, respectively, for the epoch plotted in Figure [T](a). 
Comparing these values with those from the later epochs 
of the HI 2 bxO model implies that the main effect of the 
unphysical inflow in the JM06 model is to increase the 
density of the hot halo gas by an order of magnitude, 
and hence the X-ray surface brightness by two orders of 
magnitude.) 

The predictions from model bx50 do not oscillate, but 
instead there is a general increase in the predicted X- 
ray temperatures and surface brightnesses from t = 85 
to t = 235 Myr—the medians increase by ~0.6 x 10® K 
and an order of magnitude, respectively, over this time 
period. The increase in brightness is mainly due to an 
order-of-magnitude increase in the path length through 
the hot gas, due to a shock being driven upward through 
the halo. The path length through the hot gas contin¬ 
ues to increase beyond t = 235 Myr, but more slowly 
than at earlier epochs, as the shock slows down. The 
r.m.s. density of the hot gas in the bx50 model decreases 
from t = 135 Myr onward, because the hot gas extends 
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Figure 1. Comparison of the H12 (a) bxO and (b) bx50 fountain model predictions with the HS13 halo observations in the temperature- 
surface-brightness plane. The observations are plotted in black, while the other colors represent different epochs of the H12 models. For 
each plotted data point, the symbol indicates the medians, while the error bars indicate the lower and upper quartiles. For the observations, 
only sight lines on which the halo temperature was free to vary are included in the temperature data (see HS13). The observed surface 
brightnesses have been latitude-corrected assuming a plane parallel halo geometry, and non-detections are included at their best-fit values. 
In the bxO plot, we also show results for one epoch of the JM06 model (the latest epoch examined by HIO). 




Figure 2. Comparison of the observed and predicted distributions of halo (a) temperatures and (b) 0.5-2.0 keV surface brightnesses. In 
each plot, the bottom panel shows the observed distribution (HS13). Similarly to Figure^ only sight lines on which the halo temperature 
was free to vary are included in the observed temperature distribution, and the observed surface brightness distribution shows the latitude- 
corrected values with non-detections included at their best-fit values. The arrow in the bottom panel of plot (b) indicates the median 
upper limit on the latitude-corrected surface brightness from the sight lines on which halo emission was not detected. In the other panels, 
the solid and dashed histograms show the predicted distributions from different epochs of the H12 bxO and bx50 models, respectively. The 
colors representing the different model epochs match those used in Figure fTI 
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to greater heights, and thus includes lower-density gas. 
After t = 235 Myr, the model bx50 X-ray predictions 
are fairly steady, although there is some variation in 
the shapes of the predicted distributions. At these later 
epochs of the bx50 model, its predictions are similar to 
those from model bxO. 

Now that we have understood the H12 X-ray predic¬ 
tions in terms of the physical properties of the hot gas 
in the model domains, we can compare these predic¬ 
tions with the HS13 halo measurements. The models 
generally underpredict the median observed halo tem¬ 
perature by ~10-20%, although the predictions and ob¬ 
servations agree within the observed sight line-to-sight 
line temperature variation. However, the models sig¬ 
nificantly underpredict the halo surface brightness—the 
difference is two orders of magnitude if we compare 
medians, and the predicted upper quartiles are an or¬ 
der of magnitude less than the observed lower quartile. 
Note that there is more than enough energy available 
in the model to power the observed X-ray emission in 
principle—the SN energy injection rate in the H12 model 
is 1.1 X 10^® erg s“^ kpc“^, whereas the observed 0.5- 
2.0 keV luminosity is ~8 x 10^® erg s“^ kpc~^ (HS13). 
However, in practice, only ~10“^ of the energy from SNe 
in the model is radiated as 0.5-2.0 keV photons from the 
halo. 

Despite the minor modifications to our method for ob¬ 
taining the model predictions (Section S]), the results 
from the JM06 model are consistent with those in HIO— 
as in HIO, we find that this model matches the observed 
halo surface brightness but overpredicts the halo temper¬ 
ature (Figure [TJa)). However, as noted in Section [3l the 
X-ray predictions from this model are unreliable, due to 
an unphysical inflow of hot gas into the model domain. 

6. DISCUSSION 

Here we discuss possible reasons for the discrepancy 
between the H12 model predictions and HS13’s observa¬ 
tions. First, we consider the impact of variations in the 
SN rate fSection 16.11) . We then consider the possibility 
that we are underestimating the emission from the halo 
material in the H12 model, either because we assume 
that the halo gas is in CIE fSection 16.2|) . or because we 
are underestimating the emission from interfaces between 
hot and cold gas, due to thermal conduction not being 
included in the hydrodynamical model (Section 16.3. 1|) 
and CX not being included in the emission model (Sec¬ 
tion [5321) • We then consider the possibility that cosmic 
rays (CRs) play a role in driving material out of the 
disk, meaning that the H12 model may underestimate 
the amount of hot material in the halo fSection |6.4l) . Fi¬ 
nally, we consider the role that a more extended halo of 
hot gas, predicted by galaxy formation models, ma y pla y 
in producing the observed X-ray emission (Section 16.51) . 

6 .1. Supernova Rate 

We first consider the impact of our chosen model 
parameters, such a s the SN rate and the g a s sur- 
face mass densit y. |D e Avillez &: Breitschwerdtl (120041) 
and iJoung et al.l (f^OOl) have each e xplored varia tions 
in the SN rate in similar models; IJoung et al.l cor¬ 
respondingly varied the gas surface mass dens i ty fo l- 
lowing the Kennicutt-Schmidt law (jKennicuttl 119981) . 


iDe Avillez fc Breitschwei^ (|2004ll found that the hot 
gas filling f r action increases somewhat with SN rate; 
IJoung et al.l (I2009D found that the turbulent pressure 
and thermal pressure track each other. Both found 
that the temperature of the hot gas increases some¬ 
what in higher SN rate models with relatively little 
chan ge in the density of the hot gas ( see Figures 2 and 
3 of Ide Avillez fc Breitschwerdtl 120041 and Figure 2 of 
IJoung et al.ll2009l) . We thus suspect that an increased 
SN rate would increase the X-ray temperature and, as a 
result, the X-ray surface brightness. A higher hot gas fill¬ 
ing fraction and density would also increase the surface 
bri ghtness, by in c reasin g the emission measure. Because 
the IJoung et al.l (I2009D models are not directly compa¬ 
rable to the H12 models (see Section El), a quantitative 
estimate of this effect would require running versions of 
the H12 models with varied SN rates. This is beyond the 
scope of this paper. 

Local variations in the SN history due to the pseudo¬ 
random SN distribution may also impact the observed 
properties. The consideration of multiple time steps in a 
single model addresses this source of uncertainty to some 
extent. However, the variations in both temperature and 
emission measure are relatively small over the course of 
the runs (Figures [T] and [2|). 

6 .2. Non-equilibrium Ionization 

We now consider the possibility that the HI 2 model 
underpredicts the observed halo X-ray emission because 
we assumed that the X-ray-emitting plasma was in CIE 
when calculating the X-ray spectral predictions (Sec¬ 
tion 0]) . In reality, the plasma in the halo may be ove¬ 
rionized (i.e., the ionization temperature exceeds the 
kinetic temperature) as a result of radiative or adia¬ 
batic cooling. This would result in recombination emis¬ 
sion from cool gas which is essentially non-emissive 
if we assume CIE i Bre itschwerdt fc Schmutzled 11994 
Ide Avillez fc Breitschwerdtll2012^ . Hence, bv assuming 
CIE, we may be underestimating the X-ray emission from 
the halo plasma in the H12 model. In addition, non¬ 
equilibrium ionization (NEI) would affect the radiative 
cooling rate, which would in turn affect the temperature 
structure of the halo in the hydrodynamical models—this 
too could affect the X-ray predictions. 

It is not possible to calculate the degree of overion¬ 
ization, and hence the amount of recombination emis¬ 
sion to include in the X-ray spectral predictions, when 
post-processing the HI 2 hydrodynamical data, as La- 
grangian temperature histories are not available. In¬ 
stead, one needs to trace self-consistently the ioniza¬ 
tion evolution of the relevant elements during the course 
of the hydrody namical simulation. While such s imula- 
tions do exist (Ide Avillez fc Breitschwerdtll2012allbll . X- 
ray spectral predictions that can be compared directly 
with HS13’s observations are not currently available. 
Note that, although an overionized recombining plasma 
produces a very different emission spectrum from the CIE 
plasma models used in HS13’s XMM-Newton analysis 
(free-bound versus line emission), future predictions from 
NEI ISM models could still be compared with HS13’s 
observational results, if such predictions are first charac¬ 
terized using the method described in Section S) 

Although detailed X-ray spectral predictions for a re¬ 
combining halo plasma are not currently available, we 
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can estimate by how much taking into account NEI would 
increase the predicted X-ray surface brightness of the 
HI2 mode l. For this calc ulation, we used the code de¬ 
scribed in iSheltonI (jI998fl to follow the ionization evo¬ 
lution of a stationary parcel of plasma initially in CIE 
cooling isobarically from 3 x 10® K0 At each step in 
the calculation, the code takes into account the non¬ 
equilibrium ion populations in the plasma when calculat¬ 
ing the radiative cooling function and the emergent X-ray 
spectrum. We find that, when the plasma has cooled to 
3 X 10® K, the 0.5-2.0 keV emission is ~3,000 times as 
bright as that from a CIE plasma at the same temper¬ 
ature. However, this overionized, recombining plasma is 
^17,000 times fainter than the original T = 3 x 10® K 
CIE plasma. In the H12 model, the emission measure 
of gas with T = (2-4) x 10® K is typically similar to 
(within a factor of ^^5) the emission measure of gas with 
r > 1 X 10® K0 This calculation therefore implies that 
the overionized cooled halo plasma would be much fainter 
than the hot (T > 10® K) halo plasma, and so taking into 
account overionization in the cooled halo plasma would 
not significantly increase the total X-ray surface bright¬ 
ness predicted by the HI 2 model. 


6.3. Emission from Interfaces 
6.3.1. Effect of Thermal Conduction 


We now explore the possibility that, because the HI2 
model does not include thermal conduction, we are un¬ 
derestimating the contribution to the emission from in¬ 
terfaces between tenuous hot gas (T > 10® K) and 
denser, cooler gas (T < 10^ K). Within such interfaces 
there exists X-ray-emissive gas that is denser, and as a 
result brighter per unit volume, than the diffuse hot gas. 
These interfaces are typically ~10-70 pc thick along the 
line of sight in the present model (note that these inter¬ 
faces are not well resolved in the halo, where the reso¬ 
lution is typically 16 or 32 pc in the hot gas). Thermal 
conduction would tend to broaden these interfaces until 
their widths are approxima tely equal to the Field length 
(jBegelman fc McKedll990ll . 


Ap = 



( 1 ) 


typically in the range ~ 1-600 pc (though for some inter¬ 
faces, Ap < 0.1 pc or >1 kpc). For approximately half of 
the interfaces in the HI 2 model, Ap exceeds the interface 
width, implying that thermal conduction would tend to 
broaden these interfaces. All other things being equal, 
increasing the width of such an interface increases the 
path length through the denser X-ray-emissive gas, and 
so including thermal conduction would be expected to 
boost the emission from interfaces. 

We investigated by how much the broadening of in¬ 
terfaces by thermal conduction could increase the X-ray 
emission by considering smoothly varying model inter¬ 
faces of width w between hot (Tp = several x 10® K) and 
cold (Tc = 10^ K) gas in pressure balance. In our inter¬ 
face model, the temperature across the interface varies 
with position x as 


T{x,w) 


Fc + Fh 
2 



( 2 ) 


where the interface center is located at a; = 0. As the 
interface is in pressure balance, the electron density is 

ne{x,w) =ne,h7Tr^^ (3) 

1 {x, W) 


where ne,h is the density in the hot gas. Temperature 
profiles for three example values of w are shown in Fig¬ 
ure |3Ka). Corresponding profiles of the 0.5-2.0 keV X- 
ray emission (normalized to rie^h), (fie/fT-e,h)^eo. 5 - 2 .o(T’)) 
where £o. 5 - 2 .o( 2 ^) is the plasma emissivity, are shown in 
Figure ISKb). As can be seen, an interface with u; > 0 
is locally brighter than a zero-width interface between 
X ~ —wl2 and x ~ -|-w/4. 

By integrating emission profiles like those shown in 
Figure [SDb) with respect to x, we can calculate the X- 
ray surface brightness of an interface described by Equa¬ 
tions © and m as a function of interface width w, 

1 r 

So.5-2.o{.w) = — J ^nl{x,w)eo.5-2.o(T[x,w])- 
nl{x,0)so.5-2.o{T[x,0])'^dx 

(4) 


where n is the number density. 


= 5.6 X 


10 ~'(T/K rA erg s~^ K~ cm~^ is the thermal con¬ 
ductivity (|Draine &: Giulia ni||1984D , and A is t he radia¬ 
tive cooling function ([Raymond fc Smithlll977L and up¬ 
dates) 0 Taking values from the midpoints of the in¬ 
terfaces in the H12 model (where the temperatures and 
densities are typically ^9 x 10^-7 x 10® K and ^2 x 10“'*- 
2 X 10“^ cm“®, respectively), we obtain Field lengths 


® While similar calculations have been carried out previousl y 
(e.g., IShapiro & Moor^ 119761 : Ide Avillez fc Breitschwerd^ 12012^ . 
the results are not presented in a form that can easily be applied 
to the 0.5—2.0 keV XMM-Newton band. 

^ Note that the flatness of the mass-weighted temperature distri¬ 
butions below log(T/K) ~ 6 indicates that there are similar quan¬ 
tities of ~3 X 10^ and x 10® K gas in the H12 model domains 
(see Figure 5 of the H12 erratum). 

In the definition of Ap in IBegelman &: McKed II1990I) . the 
second term in the denominator of Equation Q is Cm = 
max(A, r/n), where F is the diffuse heating rate. However, at 
the temperatures in the interfaces considered here, F = 0 (H12), 
and so Cm = 


(for a sight line looking perpendicular to the interface). 
Note that, in the above expression, we subtract off the 
emission from a zero-width interface, so S'o, 5 _ 2 .o('f'^) is 
the increase in surface brightness due to increasing an 
interface’s width from zero to w. Note also that, since 
the integrand in Equation (|3]) is a function of x/w (see 
Equations ([2]) and ([3])), Sq. 5 - 2.0 fx w. The surface bright¬ 
nesses obtained from Equation (j4|) are shown by the black 
curves in Figure (SKc), for three different values of Th. 

While increasing the widths of the interfaces would in¬ 
crease their surface brightness, in practice the increase in 
brightness cannot account for the discrepancy between 
the H12 predictions and the HS13 observations. This 
is because the emission from the diffuse hot gas tends 
to dominate over that from the interfaces between the 
hot gas and cooler gas, as we now demonstrate. The 
red curves in Figure [3] shows the surface brightnesses of 
uniform hot gas as functions of path length through the 
hot gas, for the same three values of Th used for the 
black curves. For Th < 5 x 10® K, the surface brightness 
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Figure 3. (a) Temperature profiles for the interfaces between hot 
and cold gas described by Equation (O, with Th = 5 X 10^ K 
and Tc = 1 X 10^ K, for three different interface widths {solid: 
w = 0 pc; dashed: w = 20 pc; dotted: w = 100 pc), (b) Profiles 
of the 0.5-2.0 keV emission for these model interfaces, where the 
plasma emissivity, £0.5—2.0? is a function of the temperature given 
by Equation I®, and the electron density, n©, is given by Equa¬ 
tion 0. The plot shows profiles for the same values of Th and Tc 
and for the same interface widths as in panel (a), (c): 0.5-2.0 keV 
surface brightnesses of these model interfaces (Equation 0; black 
curves) and of regions of uniform hot gas of temperature Tjj {red 
curves), as functions of interface width, w, and of path length 
through the hot gas, respectively. In each case, results are shown 
for Th = 3 X 10® {solid), 5 X 10® {dashed), and 10 X 10® K {dotted), 
assuming He h = 10“® cm“®. 

of an interface of a given width is less than the surface 
brightness of uniform hot gas of the same extent. As the 
regions of diffuse hot gas will likely be larger than the 
interfaces at their edges, the emission from the diffuse 
hot gas will tend to dominate. For example, consider a 
region of diffuse 5 x 10® K gas 500 pc in extent with zero- 
width interfaces at its edges. If thermal conduction were 
to increase the widths of those interfaces from zero to 


100 pc, the total X-ray surface brightness would increase 
by only ~20%. Therefore, increasing the widths of the 
interfaces in the model (by including thermal conduc¬ 
tion) would not counteract the two-order-of-magnitude 
difference in brightness between the HI 2 predictions and 
the HS13 XMM-Newton observations. 


6.3.2. Charge Exchange 

The above discussion considered only emission result¬ 
ing from collisional excitation of the gas in an interface. 
However, CX reactions between ions from the hot side of 
an interface and neutral H and He atoms from the cold 
side could also contribute to the emission. To estimate 
the i mportance of CX emission, we used Equation (1) 
from lLallemend (I2004f) . This gives the path length, Lex, 
through a hot gas for which the thermal emission from 
the hot gas is equal in brightness to the CX emission 
from the two interfaces at either end of the hot gas. As¬ 
suming the interfaces are observed at normal incidence, 
this path length is 


Lex = 0.06 


earicVioo 

Xnl 


pc. 


(5) 


where e is the ratio between the CX probability and the 
collisional ionization probability in the hot gas, a is the 
ratio of the global emis s ivity C X cross-section, E, to that 
assumed by iLallementl (l2004f) (E = 6 x 10“^® keV cm®, 
appropriate for solar wind CX emission in the 0.1- 
0.5 keV band), y is the ratio of the hot gas emissivity to 
that assumed bv ILallementl (5.8 x 10“^^ keV cm® s“®), 
rtc and rig are the number densities of the cold and hot 
gas, respectively, and Vioo is the relative speed of the 
ions and the neutrals in units of 100 km s“®. 

To estimate Lex, we assumed that the hot gas has a 
temperature of 3 x 10® K, and used e = 0.2 (from the 
V = 100 km curve in Figure 1 of lLallementJl200J l and 
X = 0.17 (the 0.5-2.0 keV emissivity of a 3x 10® K plasma 
is 9.8 X 10“®® keV cm® s~®: [Raymond fc SmithlflOT^ and 
updates). In the absence of suitable CX emission data 
for the XMM-Newton band, we assumed a iH We 
used a typical hot gas density of 5 x 10“® cm“® (Table (TJ 
column 5), and we assumed that the hot gas is in pressure 
equilibrium with cold gas with temperature 3 x 10^ K 
(i.e., ric = lOOrie). Finally, we assumed that the relative 
motion of the ions and neutrals is dominated by thermal 
motion, and so used Vioo 1 (cf. the mean speed of 
oxygen ions in a 3 x 10® K plasma is 60 km s“®). 

Using the above values in Equation ([5]), we find Lex 
140 kpc, i.e., a region of hot (T = 3 x 10® K, rig = 
5 X 10“® cm“®) gas in pressure equilibrium with cooler 
(3 X10"® K) gas would have to be ^140 kpc in extent in or¬ 
der for its 0.5-2.0 keV thermal emission to be as bright as 


Although we are considering a higher-energy band than 
ILallemen^ II2004I ') (0.5-2.0 versus 0.1-0.5 keV), and hence CX emis¬ 
sion from a different set of lines, we are assuming here that the 
sum of the abundances of the relevant ions and the typical CX 
cross-secti on and line yield are similar to the values that yielded 
ILallemen^ ’s assumed value of S. The line energies will of course 
be higher in the band that we are considering, which would tend 
to increase S. However, as the CX emission in the XMM-Newton 
band is likely dominated by oxygen Ka emission near 0.6 keV, the 
typical line energies in the two bands will be within a factor of a 
few of each other, and so a 1 should be a reasonable assumption. 
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the CX emission from the interfaces bounding the gas. In 
contrast, the path lengths through the hot gas in the H12 
model are typically a few kiloparsecs (Table[Tl column 7), 
implying that CX emission may be up to two orders of 
magnitude brighter than the thermal emission from the 
hot gas. It is therefore possible that CX emission could 
account for much of the shortfall between the current 
predictions from the H12 model and the observed halo 
surface brightness. However, from this simple estimate 
we cannot definitively conclude that most of the observed 
halo emission is due to CX—more detailed spectral cal¬ 
culations are needed to determine how much CX emis¬ 
sion the H12 model produces. These calculations would 
have to be carried out for each hot-cold interface in the 
model individually, taking into account the temperature 
of the hot gas (which affects the populations of the ions 
undergoing CX reactions) and the densities of the hot 
and cold gas (which affect the overall brightness of the 
CX emission), and using CX cross-section and line yield 
data suitable for emission in the XMM-Newton band. 
Such calculations are beyond the scope of this paper. 

If CX emission is indeed a major contributor to the 
observed halo X-ray emission, this would mean that the 
emission measure of the hot halo gas is smaller than pre¬ 
viously thought (e.g., ^(0.4-7) x 10“^ cm“® pc; HS13). 
This would have important implications for the results of 
joint emission-absorption analyses of the halo, in which 
emission measurements are combined with ion column 
density measurements to infer the density and extent of 
the halo. In such an analysis, the extent of the hot halo 
scales as N'^/S, where N and £ are the column density 
and emission measure of the hot gas, respectively. The 
density scales as £/N, and so for a spherical halo, the 
gas mass scales as £/N x {N^/£)^ = N^/£^. Therefore, 
if the presence of CX means that the halo emission mea¬ 
sure is overestimated, the extent and mass of the hot halo 
inferred from joint emission-absorption analyses will be 
underestimated. 

6.4. Cosmic Ray Driving 

In the HI2 model, material is driven from the disk into 
the halo solely by the thermal pressure of SN-heated gas. 
However, CRs may also play an important role in drivin g 
outflows frorn galact ic disks (jBreitschwerdt et al1ll991h . 
lEverett et al.l (1200811 showed that a CR-driven galactic 
wind (modeled in one dimension) provided a better fit 
to the diffuse 3/4 keV emission observed toward the in¬ 
ner Gal axy (—30° < I < 30°) than a static polytropic 
model. iSalem fc BrvanI (1201411 . meanwhile, used three- 
dimensional AMR simulations to study CR-driven out¬ 
flows. They showed that CR driving led to significant 
baryonic mass loss from the disk of their model galaxy, 
in contrast to a model without CR dr iving, in which 
there was no such mass loss. In addition. ISalem fc BrvanI 
(|2014ll showed that including CR diffusion (as opposed 
to just having the CRs advect along with the gas flow) 
was important for driving the outflow from the disk. 

In the context of the present study, including CR driv¬ 
ing would be expected to result in more material be¬ 
ing transported from the disk into the halo than in the 
H12 model, thus potentially increasing the halo’s X-ray 
surface brightness. However, the X-ray emission from 
such a CR-driv en outflow a l so dep ends on its tempera¬ 
ture structure. iBooth et al.l (|2013r) found that CR driv¬ 


ing results in cooler outflows than pure thermal-pressure 
driving, but they chose a feedback i mplementation equiv - 
alent to the “energy only” runs in lAgertz et al.l (l2013f) , 
which minimizes or eliminates hot gas pro duction by SN 
explosions in dense gas (see Figure 6 in lAgertz et al.f) . 
As a result, their prediction is only a lower limit on the 
true temperature. A model similar to the H12 model 
that incorporates CRs is currently under development 
(P. Girichidis et al. 2014, in preparation). The X-ray 
predictions from this new model will help determine the 
role of CR driving in supplying the hot halo gas observed 
in emission. 

6.5. Role of an Extended Galactie Halo 

Finally, we consider how an extended halo of hot gas 
(>100 kpc in extent) might affect the H12 model predic¬ 
tions. The emission from the H12 model comes mostly 
from within a few kiloparsecs of the Galactic midplane. 
This is in part due to the low densities far above the 
disk in the model (typically <3 x 10“^ cm“^ above 
10 kpc). In contrast, there is indirect evidence (from 
the lack of gas in satellite galaxies and the confinement 
of high-velocit y clouds) for hig her density halo gas far 
from the disk (iFang et al.ll2013t and references therein). 
For example, a model of an extended non-isothermal 
halo in hydrostatic equilibrium with the Galaxy’s dark 
matter (jMaller fc Bulloc3l2004f) , which is consistent with 
the observed X-ray emission and pulsar dispersion mea- 
sure data, and wit h the aforementioned indirect evidence 
(IFang et al.l 120131) . has a density exceeding 10~^ cm~^ 
out t o '^100 kpc (mb model in Figure 1 of IFang et al.l 
I2013II . Such extended hot halos a re also predicted by 
disk galaxy formation models ('e.g.. iCrain et al.ll2010ll . 

If the Milky Way’s extended halo consists of low- 
metallicity material accreted from the intergalactic 
medium, then in itself it would not be X-ray bright. How¬ 
ever, from smoothed particle hydrodynami c s (SPH ) sim¬ 
ulations of galaxy formation, ICrain et al.l (|2013f) found 
that the X-ray emission from their model galactic halos 
was produced by metals that were transported out of the 
ISM being collisionally excited by hot electrons in low- 
metallicity accreted gas. Hence, a hot, low-metallicity 
halo of accreted material could boost the X-ray emis¬ 
sion from the fountains in the HI 2 model, by increasing 
the population of electrons available to excite the ions 
in the fountains. In addition, if there have been pre¬ 
vious episodes of starburst activity in the Milky Way, 
these could have enriched the extended halo with met¬ 
als, potentially making the extended halo intrinsically 
X-ray emissive. The arbitrary inflow at the boundaries 
of the JM06 model in fact raised the high-altitude densi¬ 
ties above 10“^ cm“^. As found in HIO, this did indeed 
lead to X-ray surface brightnesses comparable to the ob¬ 
served values. If the extended halo is indeed intrinsically 
X-ray emissive, its emission would have to be added to 
that predicted by the HI 2 fountain model to get the 
total predicted halo emission. Note, however, that the 
observed halo emission is p atchy, exhibiting lar ge sight 
line-to-sight line variation (lYoshino et al.l 120091 : HS13). 
An extended halo model may have difficulty explaining 
this patchiness. 

Predictions from hydrodynamical models of galaxy for¬ 
mation are needed to test the role played by an extended 
halo in producing the X-ray emission observed from the 
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Milky Way’s halo. We plan to examine such predictions 
in a subsequent paper. 

7. SUMMARY AND CONCLUSIONS 

We have compared the X-ray emission predictions of 
a magnetohydrodynamical model of the SN-driven ISM 
(H12) with XMM-Newton measurements of the Galac¬ 
tic halo’s emission (HS13). This model significantly un¬ 
derpredicts the halo’s X-ray surface brightness (by two 
orders of magnitude, when we compare the medians of 
the predicted and observed values; Section [5]). Including 
an interstellar magnetic field does not significantly affect 
these X-ray predictions. 

We explored possible reasons for the discrepancy be¬ 
tween the H12 model predictions and HSlS’s XMM- 
Newton observations. Assuming CIE may in principle 
underestimate the emission from the H12 model halo, 
but in practice this is unlikely to have a significant effect 
(Section 15^ . We also found that the discrepancy could 
not be explained by the emission from interfaces in the 
HI2 model being underestimated due to a l ack of ther¬ 
mal conduction in the model (Section 16. 3. II) . However, 
CX emission from such interfaces (not included in the 
present emission model) could greatly increase the pre¬ 
dicted X-ray surface brightness, though detailed s pectra l 
calculations are needed to confirm this (Section l6.3.2p . 
(If CX emission is a major contributor to the observed 
halo emission, then the hot gas emission measure is less 
than previously thought, with the consequence that the 
path length and mass of the hot gas calculated from al¬ 
gebraic combinations of the emission measure and ion 
column density would be revised upwards.) In addition, 
CR driving of a wind could increase the amount of X- 
ray-emissive material in the halo (Section 16.41) . and an 
extended hot halo of accreted material, while not intrin¬ 
sically X-ray bright, may supply hot electrons that could 
increase the predicted X-ray emission from galactic foun¬ 
tains (Sectioning. 

In conclusion, the faintness of the HI 2 model relative 
to the observed surface brightness implies that thermal 
emission from classical galactic fountains is not a major 
source of the halo’s X-ray emission. This is in contrast 
to the conclusion of HIO, which was based on the JM06 
model (the X-ray predictions from which are now known 
to be incorrect). Our results indicate that additional 
physical processes need to be included in halo models. 
Two plausible possibilities are the effects of CR driving 
on the fountain, and extended hot halos from the galaxy 
formation process. In addition, CX may be an impor¬ 
tant contributor to the observed emission. Suitable X- 
ray predictions from CR-driven ISM models and galaxy 
formation models are needed to test the roles of galactic 
fountains and of accreted extragalactic material in ex¬ 
plaining the observed X-ray emission from the Galactic 
halo. 
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